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The Polya number characterizes the recurrence of a random walk. We apply the generalization 
of this concept to quantum walks [M. Stefanak, I. Jex and T. Kiss, Phys. Rev. Lett. 100, 020501 
(2008)] which is based on a specific measurement scheme. The Polya number of a quantum walk 
depends in general on the choice of the coin and the initial coin state, in contrast to classical random 
walks where the lattice dimension uniquely determines it. We analyze several examples to depict the 
variety of possible recurrence properties. First, we show that for the class of quantum walks driven 
by tensor-product coins, the Polya number is independent of the initial conditions and the actual 
coin operators, thus resembling the property of the classical walks. We provide an estimation of 
the Polya number for this class of quantum walks. Second, we examine the 2-D Grover walk, which 
exhibits localisation and thus is recurrent, except for a particular initial state for which the walk is 
transient. We generalize the Grover walk to show that one can construct in arbitrary dimensions a 
quantum walk which is recurrent. This is in great contrast with classical walks which are recurrent 
only for the dimensions d = 1, 2. Finally, we analyze the recurrence of the 2-D Fourier walk. This 
quantum walk is recurrent except for a two-dimensional subspace of the initial states. We provide 
an estimation of the Polya number in its dependence on the initial state. 
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I. INTRODUCTION 

Random walks (RWs) present a very useful tool in 
many branches of science yj]. Among others, the RW 
is one of the cornerstones of theoretical computer science 
[3, Q. Indeed, it can be employed for algorithmic pur- 
poses solving problems such as graph connectivity 
n-SAT or approximating the permanent of a matrix 
0. 

Quantum random walks (QWs) as a generalization of 
classical random walks to quantum systems have been 
proposed by Aharonov, Davidovich and Zagury The 
unitary time evolution can be considered cither dis- 
crete as introduced by Meyer 0, [§] and Watrous 
leading to coined QWs or continuous as introduced by 
Farhi and Gutman (Til . [T^ . Scattering quantum walks 
[isl fl3 . m, nil were proposed by Hillery, Bergou and 
Feldman as a natural generalization of coined QWs based 
on an interferometric analogy. The connection between 
the coined QW and the continuous time QW has been 
established recently [17 1. 

The coined QW is well suited as an algorithmic tool 
[Tsl . [l9| . Several algorithms based on coined QWs show- 
ing speed up over classical algorithms have been pro- 
posed 0, dll, m, [2^ . Various properties of coined 
QWs have been analyzed, e.g. the effects of the coin 
and the initial state [25l. |26|. Wf\ , absorbing barriers [1^ , 
the hitting times 2^, 30l l3l| or the effect of decoherence 
[U [S^l- Hitting times for continuous QWs related to 
the quantum Zeno effect were considered in [S^l- Great 
attention has also been paid to the asymptotics of QWs 
[3^. Issl. IsgI. IstI. Isst . In particular, localisation was found 



in 2-D QWs [25|, l39|, lid] and in 1-D for a generalized QW 
[41I [i^ l- Several experimental schemes have been pro- 
posed to realize coined QWs including cavity QED 43 1, 
linear optics [44, 45], optical lattices 46, 47] to BEG [48| . 
By now, quantum walks form a well established part of 
quantum information theory (49| . 

An interesting question for classical random walks is 
whether the particle eventually returns to the starting 
point of the random walk. The recurrence probability 
is known as the Polya number, after G. Polya who first 
discussed thisproperty in classical RWs on infinite lat- 
tices in 1921 [50|. Polya pointed out the fundamental dif- 
ference between walks in different dimensions. In three 
or higher dimensions the recurrence has a finite, non- 
unit probability depending exclusively on the dimension, 
whereas for walks in one or two dimensions the Polya 
number equals one. As a consequence, in three and 
higher dimensions the particle has a non-zero probability 
of escape [5l|, HI] . Recurrence in classical RWs is closely 
related to first passage times [s^. Recurrence in semi- 
infinite and finite QWs on lattices has been calculated in 

[l,!!!. 

In a recent letter [HsI] we have defined the Polya number 
for quantum walks on a d dimensional lattice by extend- 
ing the concept of recurrence. In this paper we calculate 
the Polya number for various coined QWs in one and two 
dimensions and construct arbitrary dimensional QWs ex- 
hibiting highly non-classical features. 

Our paper is organized as follows: In Section |TT] we 
review the concept of recurrence and the Polya number 
of random walks and its extension to quantum walks as 
defined in ^55']. Both classical and quantum definitions 
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lead to a similar criterion for recurrence determined by 
the asymptotic scaling of the return probability po(i)i 
we show in Appendix [Al 

We dedicate Section HTll to the study of the asymptotic 
behaviour of the return probability po(t). For this pur- 
pose we apply the Fourier transformation and the method 
of stationary phase. In particular, we demonstrate that 
the asymptotic behaviour of the return probability po (t) 
is influenced by three factors: the topology of the walk, 
the choice of the coin operator and the initial coin state. 
Consequently, the nature of the QW for a fixed dimen- 
sion can change from recurrent to transient and the ac- 
tual value of the Polya number varies. This is in great 
contrast to classical random walks where the recurrence 
is uniquely determined by the dimension of the walk. 

We use the results derived in Section Hill to determine 
the recurrence properties of several types of QWs. In 
Section IIV Al we treat unbiased 1-D QWs and show that 
all of them are recurrent independent of the coin operator 
or the initial coin state. We then generalize 1-D QWs in 
Section [IV Bl to d dimensions by considering an indepen- 
dent coin for each spatial dimension. We find that for this 
class of QWs the Po{t) is independent of the initial coin 
state and the actual form of the coin operator. Hence, 
a unique Polya number can be assigned to this class of 
QWs for each dimension d. In contrast with the classical 
RWs this class of QWs is recurrent only for d = 1. 

In Section IV Al we analyze the recurrence of the 2-D 
Grover walk. This QW exhibits localisation [4^ and 
therefore is recurrent. However, for a particular initial 
state localisation disappears and the QW is transient. 
We find an approximation of the Polya number for this 
particular initial state. In Section lVBI we employ the 2-D 
Grover walk to construct for arbitrary dimension d a QW 
which is recurrent. This is in great contrast with the clas- 
sical RWs, which are recurrent only for the dimensions 
d= 1,2. 

Finally, in Section lVTl we analyze the 2-D Fourier walk. 
This QW is recurrent except for a two-parameter family 
of initial states for which it is transient. For the latter 
case we find an approximation of the Polya number de- 
pending on the parameters of the initial state. 

We conclude by presenting an outlook in Section IVlII 

II. RECURRENCE AND THE POLYA NUMBER 
OF RANDOM WALKS 

Random walks are classically defined as the probabilis- 
tic evolution of the position of a point-like particle on a 
discrete graph. Starting the walker from a well-defined 
graph point (the origin) one can ask about the probabil- 
ity that the walker returns there at least once during the 
time evolution. The event that the walker is not present 
at the origin at any time instant is just the complement 
of the event corresponding to recurrence. The probabil- 
ity of the latter is called the Polya number. Classical 
random walks are said to be recurrent or transient de- 



pending on whether their Polya number equals to one, 
or is less than one, respectively. 

The Polya number of a classical random walk can be 
defined in the following way 56]. Let qait) be the prob- 
ability that the walker returns to the origin for the first 
time after t steps. Since these events are mutually exclu- 
sive we can add up their probabilities and the series 

oo 

t=i 

is the probability that at least once the particle has re- 
turned to the origin, i.e. the Polya number. However, 
the definition ([T]) is not very practical for determining 
the recurrence nature of a random walk. Nevertheless, 
we can express the Polya number in terms of the proba- 
bility Po (t) that the particle is at the origin at any given 
time instant t. Indeed, it is straightforward to show that 



+00 ^ ^ 

E po{t) 

t=0 

Hence, the recurrence behaviour of a RW is determined 
solely by the infinite sum 

00 

5^EpoW. (3) 

t=o 

We find that P equals unity if and only if the series S 
diverges [HBl- In such a case the walk is recurrent. On the 
other hand, if the series S converges, the Polya number 
P is strictly less than unity and the walk is transient. 
We will restrict the considered graphs to d dimensional 
uniform infinite square lattices and the considered walks 
to balanced ones, i.e. the probability to take one step is 
equal in each of the possible directions along the lattice. 
Polya proved in that in one and two dimensions, such 
RWs are recurrent while for d > 2 they are transient. 

Discrete time coined quantum walks are generaliza- 
tions of classical random walks. Here the dynamics is 
a unitary evolution on a composite Hilbert space consist- 
ing of the external positions on the graph and the internal 
state determining the direction of the following step. For 
quantum walks one would like to keep the same definition 
of the Polya number, i.e. the probability of returning to 
the origin at least once during the time evolution. How- 
ever, to be able to talk about the position of a particle 
in quantum mechanics one must specify when and which 
type of measurement is performed. We have given a def- 
inition of the Polya number for QWs [SSj by specifying 
the following measurement scheme. From an ensemble of 
identically prepared QW systems take one, let it evolve 
for one step, measure the position and discard this sys- 
tem; take a second, identically prepared system, let it 
evolve for two steps, measure the position and then dis- 
card it; continue until it is found at the origin. In the 
t-th trial we do not find the particle at the origin with 
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the probability 1 — po (t) . Since the individual trials are 
independent the product 

n 

Pn = l[{l-P0{t)) (4) 

t=l 

gives the probability that we have not found the particle 
at the origin in the first n trials. In the complementary 
event which occurs with the probability 

n 

Pn = l-l[(^-P0{t)) (5) 
t=l 

the particle was found at least once at the origin, by 
which we understand the recurrence. Hence, we define 
the Polya number of a QW by extending n to infinity 

P = l-l[{l-po{t)). (6) 

t=l 

As we show in Appendix \^ Eq. ([6]) leads to the same 
criterion for recurrence in terms of po{t) — the infinite 
product in Eq. ^ vanishes if and only if the series S 
diverges [13] ■ In such a case the Polya number of a QW 
is unity and we call such QWs recurrent. If the series S 
converges, then the product in Eq. ^ does not vanish 
and the Polya number of a QW is less than one. In accor- 
dance with the classical terminology we call such QWs 
transient. We use the formula in Eq. ([5]) to estimate the 
Polya number of a QW. 

To conclude this section, the recurrence nature of 
quantum random walks is determined by the decay of 
the return probability po(t), in the same way as for clas- 
sical random walks. 



III. ASYMPTOTICS OF QUANTUM WALKS 

Before we turn to the recurrence of QWs we analyze 
the asymptotics of the return probability po (t) . Our ap- 
proach is based on the Fourier transformation and the 
method of stationary phase. We consider random walks 
where the walker has to leave its actual position at each 
step. Hence, po(t) = for odd times and it is sufficient 
to consider only even times 2t. For simplicity we omit 
the factor of 2 from now on. 

A. Description of quantum walks on Z'' 

We consider quantum walks on an infinite d dimen- 
sional lattice Z"^. The Hilbert space of the quantum walk 
can be written as a tensor product 

n^rLp®nc (7) 

of the position space 

np=£'^{Z'^) (8) 



and the coin space Tic- The position space is spanned 
by the vectors |m) corresponding to the walker being at 
the lattice point m, i.e. 

Tip = Span{|m)| m {mi, . . . , Wd} e Z''} . (9) 

The coin space He is determined by the topology of the 
walk. In particular, its dimension c is given by the num- 
ber of possible displacements in a single step. We denote 
the displacements by vectors 

e.eZ"^, i = l,...,c. (10) 

Hence, the walker can move from m to any of the points 
m-|-ei,i = l,...,cina single step . We restrict ourselves 
to unbiased walks where the displacements satisfy the 
condition 

c 

^e, = 0. (11) 

1=1 

We define an orthonormal basis in the coin space by as- 
signing to every displacement e.; the basis vector je^), i.e. 

Tic- Span {|e,)|i = l,...,c}. (12) 

A single step of the QW is given by 

U^S{Ip(E)C). (13) 

Here Ip denotes the unit operator acting on the position 
space Tip. The coin flip operator C is applied on the 
coin state before the displacement S itself. The coin flip 
C can be in general an arbitrary unitary operator acting 
on the coin space Tic- We restrict ourselves to unbiased 
walks for which the coin C meets the requirement 

|C.,|^|(e,|C|e,)H^, (14) 

i.e. all matrix elements of C must have the same absolute 
value. Such matrices are closely related to the Hadamard 
matrices fssj . 

The displacement itself is represented by the step op- 
erator S 

S' = ^|nH-e,>(m|(8)|e,)(e,|, (15) 

which moves the walker from the site m to m + if the 
state of the coin is \ei). 

Let the initial state of the walker be 

1^(0)) =^^,(m,0)|m)® |e,). (16) 

m, t 

Here ?Ai(m, 0) is the probability amplitude of finding the 
walker at time i = at the position m in the coin state 
\ei). The state of the walker after t steps is given by suc- 
cessive application of the time evolution operator given 
by Eq. (|13p on the initial state 

|^(t)) =^^,(m,t)|m)® |e,) = (7*|7^(0)). (17) 
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The probability of finding the walker at the position m 
at time t is given by the summation over the coin state, 
i.e. 



pim,t) EE ^|(m|(e,|^(t))p=^|^,(m,i)|2 

= ||^(m,t)|p. (18) 
Here we have introduced c component vectors 

^im,t) = {iMm,t),ib2im,t), . . . ,Mm,t)f (19) 

of probability amplitudes. We rewrite the time evolution 
equation pT|) for the state vector \tp{t)) into a set of 
difference equations 



ip{m,t) = ^C/V'(m-e/,<- 1) 



(20) 



for probability amplitudes 4'{^^t)- Here the matrices 
Ci have all entries equal to zero except for the l-th row 
which follows from the coin-flip operator C, i.e. 



(e^ \Ci\ej) = Sii{ei \C\ ej) 



B. Solution via Fourier Transformation 



(21) 



The QWs we consider are translationally invariant 
which manifests itself in the fact that the matrices C/ 
on the right-hand side of Eq. (PP]) are independent of 
m. Hence, the time evolution equations (|^n|) simplify 
considerably with the help of the Fourier transformation 



(22) 



The Fourier transformation defined by Eq. is an 

isometry between P{'L'^) and L^(K'*) where K = (— tt, tt] 
can be thought of as the phase of a unit circle in R^. 

The time evolution in the Fourier picture turns into a 
single difference equation 



V'(k,t) = c/(k)vXk,t-i). 



(23) 



Here we have introduced the time evolution operator in 
the Fourier picture 



t/(k) = Z3(k)C 
D{\i) = diag(e-*'=i''. 



(24) 



We find that ?7(k) is determined both by the coin C 
and the topology of the QW through the diagonal matrix 
Z?(k) containing the displacements e^. 

We solve the difference equation ([^5]) by formally di- 
agonalising the matrix t/(k). Since it is a unitary matrix 
its eigenvalues can be written in the form 



Aj(k) — exp (zcL'j(k)). 



(25) 



We denote the corresponding eigenvectors as Wj(k). Us- 
ing this notation the state of the walker in the Fourier 
picture at time t reads 

V;(k,i) ^ ^e-^-(k)* (v;(k,0),t;,(k)) T;,(k), (26) 



where ( , ) denotes the scalar product in the c dimen- 
sional space. Finally, we perform the inverse Fourier 
transformation and find the exact expression for the 
probability amplitudes 



i/'(m,t) 



rfk 

(27r)' 



V'(k,t) 



(27) 



in the position representation. 

We are interested in the recurrence nature of QWs. As 
we have shown in Section|II]the recurrence of a QW is de- 
termined by the asymptotic behaviour of the probability 
that the walker returns to the origin at time t 



(28) 



Hence, we set m = in Eq. (P7)) . Moreover, in analogy 
with the classical problem of Polya we restrict ourselves 
to QWs which start at origin. Hence, the initial condition 
reads 

^{m,0) = Sm^o^j, V = V'(O,0) (29) 

and its Fourier transformation ipO^j 0) entering Eq. (|26p 
is identical to the initial state of the coin 

V^(k,0)=^, (30) 

which is a c-component vector. We note that due to the 
Kronecker delta in Eq. (|29p the Fourier transformation 
'0(k, 0) is a constant vector. 

Using the above assumptions we find the exact expres- 
sion for the return probability 



Po{t) = 



i=i 



(31) 



where Ij (t) are given by the integrals 



/,(k) = (V',«,(k)) v,{k). 



(32) 



C. Asymptotics of pQ{t) via the method of 
stationary phase 

Let us now discuss how the additional freedom we have 
at hand for QWs influences the asymptotics of the re- 
turn probability po{t). For simplicity we suppose that 
the asymptotic behaviour of pQ{t) arises from \Ij{t)\'^. 
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We suppose that the functions ^jjCk) and /j(k) entering 
Ij (t) are smooth. According to the method of stationary 
phase [5§] the major contribution to the integral Ij{t) 
comes from the saddle points k° of the phase Wj(k), i.e. 
by the points where the gradient of the phase vanishes 



k=kO 



0. 



(33) 



The asymptotic behaviour of Ij (t) is then determined by 
the saddle point with the greatest degeneracy given by 
the dimension of the kernel of the Hessian matrix 



(34) 



evaluated at the saddle point. The function /j(k) en- 
tering the integral Ij{t) determines only the pre- factor 
of the asymptotic behaviour. We now discuss how does 
the existence, configuration and number of saddle points 
affect the asymptotic behaviour of Ij{t). As a rule of 
thumb, the decay of the return probability po(^) can slow 
down with the increase in the number of saddle points. 
We briefly summarize the cases of no saddle points, fi- 
nite number of saddle points and a continuum of saddle 
points. 

(i) No saddle points 

If ujj (k) has no saddle points then Ij [t) decays faster 
than any inverse polynomial in t. Consequently, the de- 
cay of the return probability po (0 would also be exponen- 
tial. However, among the examples we have considered 
such a situation was not found. 

(?i) Finite number of saddle points 

If ujj (k) has a finite number of non-degenerate saddle 
points, i.e. the determinant of the Hessian matrix H 
is non-zero for all saddle points, and the function /j(k) 
does not vanish at the saddle points then the contribu- 
tion from all saddle points to the integral Ij{t) is of the 
order and po(i) ~ t "^- Though the contributions 

from the distinct saddle points might have different rela- 
tive phases and can interfere destructively we have never 
encountered a complete cancelation of all contributions. 

(m) Continuum of saddle points 

If ojj (k) has a continuum of saddle points then the di- 
mension of the continuum determines the decay of the 
integral Ij{t). The case of 2-D integrals with curves of 
stationary points are treated in [s^ . It is shown that the 
contribution from the continuum of stationary points to 
the integral Ij{t) is of the order t~^/'^. This is greater 
than the contribution arising from a discrete saddle point 
which is of the order . Hence, the continuum of sad- 
dle points has effectively slowed-down the decay of the 
integral Ij{t). Consequently, the leading order term of 
the return probability is pa{i) ^ . Similar results 
can be expected for higher dimensional QWs where the 
phase ijJji)^) has a continuum of saddle points. A spe- 
cial case for a continuum of saddle points is when ujj (k) 
does not depend on n variables but it has a finite num- 
ber of saddle points with respect to the remaining d — n 
variables. Indeed, such an Wj(k) has obviously a zero 



derivative with respect to those n variables. If, in addi- 
tion, Ij (t) factorizes into the product of time- independent 
and time-dependent integrals over n and d — n vari- 
ables and the time-independent integral does not vanish 
then Ij{t) behaves asymptotically like Hence, 
the asymptotic behaviour of the return probability is 
Po{t) ^ t~'-'^~"\ In the extreme case when the phase 
ti;j(k) does not depend on k at all we can extract the 
time dependence out of the integral Ij (t) . If the remain- 
ing time independent integral does not vanish then pq (t) 
is a non-zero constant. 

In the above discussion we have assumed that the func- 
tion fj (k) is non- vanishing for k values corresponding to 
the saddle points. However, the initial state ip can be 
orthogonal to the eigenvector Vj(k) for k = k° corre- 
sponding to the saddle point. In such a case the function 
/j(k) vanishes for k = k" and the saddle point k° does 
not contribute to the integral Ij{t). Consequently, the 
decay of poit) can speed up. 

In Section |TT] we have shown that the recurrence nature 
of the QWs is determined by convergence or divergence of 
the sum ([3]) which in turn depends on the speed of decay 
of the return probability pQ{t). Hence, for QWs we might 
change the recurrence behaviour and the actual value of 
the Polya number by altering the initial state V': coin 
flip C and the topology of the walk determined by the 
displacements e^. 

In the following sections we make use of the above de- 
rived results and determine the recurrence behaviour and 
the Polya number of several types of QWs. We concen- 
trate on the effect of the coin operators and the initial 
states. For this purpose we fix the topology of the walks. 
We consider QWs where the displacements have all 
entries equal to ±1 



ei = (1,...,!)' 



,e2d = (-1, . 



(35) 



In such a case the coin space has the dimension = 2'^ 
and the diagonal matrix D(k) can be written as a tensor 
product 



D{k)^D{h) 



■ ®D{kd) 
- diag(e~*'^' 



(36) 
■*'='). This 



of 2 X 2 diagonal matrices D{kj) — diag^e '"^ , e 
fact greatly simplifies the diagonalisation of the time evo- 
lution operator in the Fourier picture U (k) . 



IV. RECURRENCE OF 1-D QWS AND QWS 
WITH TENSOR-PRODUCT COINS 



We begin this section with the analysis of unbiased 
1-D QWs. We find that all unbiased 1-D QWs are recur- 
rent independently of the initial coin state and the actual 
form of the coin operator. We then generalize unbiased 
1-D QWs to d dimensions by considering coins which can 
be written as a tensor products of d 2 x 2 matrices, i.e. 
we consider independent coin for each spatial dimension. 



6 



This class of d dimensional QWs maintains some prop- 
erties of the 1-D QWs. In particular, the asymptotic 
behaviour of the probability po(^) is independent of the 
initial coin state and the actual form of the coin opera- 
tor. Hence, a unique Polya number can be assigned to 
this class of QWs for each dimension d. In contrast with 
the classical RWs they are recurrent only for d = 1. 



A. Unbiased QWs in one dimension 



It follows that also the time evolution operator in the 
Fourier picture described by Eq. ([M]) has the form of the 
tensor product 

U^^^ (k, a, /3) - U{h , ai, ® . . . ® C7(fcd, a^, Pd) (42) 



of d 1-D time evolution operators given by Eq. ([38]) with 
different parameters ki,ai^Pi. Hence, the phases of the 
eigenvalues of the matrix [/'^'^^(k, a,/?) have the form of 
the sum 



Let us start with the analysis of the recurrence be- 
haviour of unbiased 1-D QWs. The general form of the 
unbiased coin for 1-D quantum walk is given by 



C(a,/3) 



V2 ( 



-1/3 



,if3 



(37) 



We find that the time evolution operator in the Fourier 
picture 



U{k,a,(3) = D{k)C{a,f3) 



(38) 



has eigenvalues e^'^iik,a-) -^^jth the phases U!i{k,a) given 
by 

n \ sin(fc-Q;) 
smuji{k,a) = p , a;2(fc, a) = tt — a;i(/e, a). 



V2 

Thus the derivatives of Ui with respect to k reads 
duji(k,a) dijj2{k,a) cos (fc — a) 



dk 



dk 



2 - sin^ (fc - a) 



(39) 
(40) 



and we find that the phases a;i(fc, a) have common non- 
degenerate saddle points k'^ — a ± tt /2. It follows that 
PQ^t) behaves asymptotically like t"^, independently of 
the coin parameters a, (3. Moreover, the asymptotic be- 
haviour is independent of the initial state. Indeed, no 
non-zero initial state if) exists which is orthogonal to both 
eigenvectors at the common saddle points fc° = a ± 7r/2. 
Hence, all unbiased 1-D quantum walks are recurrent, i.e. 
the Polya number equals one, independently of the ini- 
tial coin state and the coin. However, none of the QWs 
from the class described by Eq. ([57|) exhibits localisa- 
tion, since for all of them the probability po {t) converges 
to zero. We note that one can achieve localisation in 1-D 
by considering generalized QWs for which the coin has 
more degrees of freedom (4l| . 



B. Higher dimensional QWs with tensor-product 
coins 



We now turn to the class of QWs with independent 
coin for each spatial dimension, i.e. the coin fiip operator 
has the form of the tensor product of d 2 x 2 matrices 

C('^)(a,/3) = C(ai,/3i)®...®C(ad,/3d). (41) 



u;j(k,a) = ^Wj,(fci,ai). 



(43) 



1=1 



Therefore we find that the asymptotic behaviour of this 
class of QWs follows directly from the asymptotics of the 
1-D QWs. Indeed, the derivative of the phase a;j(k, a) 
with respect to ki reads 



9a;j(k, a) _ dujj,{ki,ai) 



dk. 



dk. 



(44) 



and so Wj(k, a) has a saddle point k° = (fcj, fc^, . . . , fc^) 
if and only if for all Z = l,...,d the point fcj* is the 
saddle point of tjj, (fc;, a;). As we have found from 
Eq. ((40|l the saddle points of LUj^ are kf — ai ± tt/2. 
Hence, all phases ujjCk, a) of the eigenvalues of the ma- 
trix L/('')(k,a,/3) have 2'* common saddle points k = 
(ai ± 7r/2, . . . , ttd ± 7r/2). It follows that the asymptotic 
behaviour of the probability pq (t) is determined by 



plf\t)^t- 



(45) 



As follows from the results for 1-D QWs the asymptotic 
behaviour given by Eq. is independent of the initial 
coin state and of the coin parameters a, (3. Compared to 
classical walks this is a quadratically faster decay of the 
probability po (t) which is due to the quadratically faster 
spreading of the probability distribution of the 1-D QWs. 

We illustrate the results for 2-D Hadamard walk driven 
by the coin which is a tensor product of two 2x2 
Hadamard matrices 



C(2)(0,0) 




(46) 



in Figure [T] Here we show the probability distribution 
and the probability po{t). The first row indicates that 
the initial state of the coin infiuences mainly the edges 
of the probability distribution. However, the probability 
Po{t) is unaffected and is exactly the same for all initial 
states. The lower plot confirms the asymptotic behaviour 
Po{t) ^ 

Since the return probability po{t) decays like i"'' we 
find that d dimensional quantum walks with tensor- 
product coins are recurrent only for dimension d = 1 and 
are transient for all higher dimensions d > 2. Moreover, 
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10- 



5 ■ 10- 



-50 




50-50 



= (i,(U), 0)^ 




arbitrary i/^ 




FIG. 1: Probability distribution of the 2-D Hadamard walk 
after 50 steps and the probability po{t) for different choices 
of the initial state. In the upper plot we choose the initial 
state ^(1,2,2,-1)"^ which leads to a symmetric probability 
distribution, whereas in the middle plot we choose the initial 
state (1,0,0,0)"^ resulting in a dominant peak of the proba- 
bility distribution in the lower-left corner of the (m, n) plane. 
However, the initial state influences the probability distribu- 
tion only near the edges. The probability po{t) is unaffected 
and is the same for all initial coin states. The lower plot con- 
firms the asymptotic behaviour po{t) ^ independent of 
the initial state. 



the whole sequence of probabihties po{t) is independent 
of the initial state and the coin C('')(a,/3). Hence, the 
Polya number for this class of QWs depends only on the 
dimension of the walk d, thus resembling the property of 
the classical walks. On the other hand, this class of QWs 
is transient for the dimension d = 2 and higher. This is a 
direct consequence of the faster decay of the probability 
at the origin which, in this case, cannot be compensated 
for by interference. 



Dimension 


Simulation 


Estimation (|49|l 


Error (%) 


2 


0.29325 


0.27325 


6.8 


3 


0.12947 


0.12841 


0.82 


4 


0.06302 


0.06296 


0.01 


5 


0.031313 


0.031309 


0.01 



TABLE I: Comparison of the Polya numbers for the class of d 
dimensional QWs with tensor product coins obtained from the 
numerical simulation and the and the analytical estimation of 
Eq. dMJ. 



Let us now estimate the value of the Polya number 
for the class of QWs with tensor-product coins and the 
dimension d > 2. As depicted in the lowest plot of Fig- 
ure [1] the return probability approaches quite rapidly its 
asymptotic 



1 



(47) 



Hence, already first few terms of the product in Eq. ([5]) 
are sufficient to estimate the value of the Polya number. 
Taking into account the first three terms of p^f^ (t) which 
are found to be 



1 



(48) 



we obtain the following approximation of the Polya num- 
ber 



pid) _ ;^ _ 



(49) 



We compare the estimation in Eq. (j49p with the nu- 
merical results obtained from the simulation of the QWs 
with 1000 steps in the Tableland find that they are in 
excellent agreement. 



V. RECURRENT QUANTUM WALKS BASED 
ON THE 2-D GROVER WALK 

We now turn to the QWs based on the 2-D Grover 
walk. This QW has been extensively studied by many 
authors [1^ [H, HO, UHl . We re-derive the properties of 
this QW using the tools developed in Section IIIII We 
find that the 2-D Grover walk exhibits localisation and 
is therefore recurrent except for a particular initial state. 
We find an estimation of the Polya number in the latter 
case. 

Employing the 2-D Grover walk we construct in arbi- 
trary dimensions a QW which is recurrent, except for a 
subspace of initial states. This is in striking contrast to 
the classical random walks which are recurrent only for 
the dimensions d = 1,2. 
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picture defined by Eq. (1241) for the Grover walk 
t/G(fci,fe) = iDih)®D{k2))G. 



We find that they are given by 
where the phase Lu{ki,k2) reads 



^±iuj{ki .ki) 



cos(Lj(fci, k2)) — — COS /ci cos k2- 



(52) 



(53) 



(54) 



V-' = i'G 




FIG. 2: Probability distribution of the Grover walk after 50 
steps for different choices of the initial state. In the upper plot 



The eigenvalues Ai,2 are constant. As a consequence the 
return probability is non-vanishing as discussed in de- 
tail in Section UlI Ci unless the initial state is orthogonal 
to the eigenvectors corresponding to Ai,2 at every point 
(/ci, ^2). By explicitly calculating the eigenvectors of the 
matrix Uciki, k2) it is straightforward to see that such a 
vector is unique and equals that in Eq. (jSTj) , in agreement 
with the result derived in [ioj . 

It is now easy to show that for the particular initial 
state given by Eq. ([5T|) the probability po [t) decays like 
t^"^ . Indeed, as the initial state of Eq. (jSTjl is orthogo- 
nal to the eigenvectors corresponding to Ai,2 the asymp- 
totic behaviour is determined by the remaining eigenval- 
ues A34(fci,fc2), or more precisely by the saddle points 
of the phase ti;(fci,A:2)- From Eq. II54I) we find that it 



we choose the initial state tps = !(!,«,«, -1)^ which leads to has only non-degenerate saddle points k^, k^ = ±n/2. 



a symmetric probability distribution with a dominant central 
spike. However, if we chose the initial state ipc according to 
Eq. (|51|l we find that the central spike vanishes and most 
of the probability is situated at the edges, as depicted in the 
lower plot. 



A. 2-D Grover walk 

We start with the 2-D Grover walk which is driven by 
the coin 



G = 




(50) 



It was identified numerically [25| and later proven ana- 
lytically [lOl that the Grover walk exhibits a localisation 
effect, i.e. the probability po(i) does not vanish but con- 
verges to a non-zero value except for a particular initial 
state 

^G = V'G(0,0,0) = i (1,-1, -1,1)^. (51) 

To illustrate this fact we present in Figure [2] the proba- 
bility distribution of the Grover walk for different choices 
of the initial coin state, namely for -05 — i(l,i,i,— 1)^ 
and for ipc given by Eq. ([5T|) . 

In order to explain the localisation we analyze the 
eigenvalues of the time evolution operator in the Fourier 



For the initial state of Eq. (|5ip the probability that the 
Grover walk returns to the origin decays like We 
conclude that the Grover walk on a 2-D lattice is recur- 
rent and its Polya number equals one for all initial states 
except the one given in Eq. (|5ip for which the walk is 
transient. 

We illustrate this result in Figure [3] where we plot the 
probability po (t) for the two different choices of the initial 
coin states, namely = 5(1, ~1)"^ and i/jg- The 
plots confirm the analytical results of the scaling of the 
probability po{t)' on the upper plot we observe that po(t) 
for the state rps oscillates around a nonzero value and 
thus has a non- vanishing limit, whereas on the lower plot 
we find that the probability po (t) for the state i/jq decays 
like t-2. 

Let us estimate the Polya number of the Grover walk 
for the initial state of Eq. (|5ip . The numerical simu- 
lations indicate that the return probability po(^) for the 
initial state "^G is the same as the return probability of 
the 2-D QW with tensor product coin. Hence, their Polya 
numbers coincide. With the help of the relation we 
can estimate the Polya number of the Grover walk with 
the initial state of Eq. (|?T|) by 

Pg(V'g) = P'-^^ ~ 0.27325. (55) 



B. Recurrent quantum walks in arbitrary 
dimensions 



The above derived results allow us to construct for an 
arbitrary dimension d a QW which is recurrent, except 
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'4> = 'ips 



0,37 



o 





o 

SSh 0.3 



FIG. 3: The probability po{t) for the Grover walk and dif- 
ferent choices of the initial coin state. First, the walk starts 
with the coin state tps ~ i, i, — 1)"^ which leads to the 
non- vanishing value of po{t), as depicted in the upper plot. 
The situation is the same for all initial coin states, up to 
the asymptotic value of po{t +oo), except for tpa ~ 
i (1, —1, —1, 1)"^, as shown in the lower plot. Here we plot 
the probability po{t) multiplied by to unravel the asymp- 
totic behavior of po{t). The plot confirms the analytic result 
of the scaling po{t) ~ t~^. 



for a subspace of initial states. Let us first consider the 
case when the dimension of the walk is even and equals 
2d. We choose the coin as a tensor product 



(56) 



of d Grover coins given by Eq. ([5(j| . As follows from Eqs. 
([M)) and the time evolution operator in the Fourier 
picture is also a tensor product 

U^G'^\k) ^UaikiM)®---® UG{k2d-i,k2d) (57) 



of the matrices Ug defined by Eq. ((52)) with different 
Fourier variables ki. Hence, the eigenvalues of [/^^''■'(k) 
are given by the product of the eigenvalues of C/q. Since 
two eigenvalues of Ug are constant as we have found in 
Eq. (|53p one half of the eigenvalues of C/^^'^^(k) are also 
independent of k. As we have discussed in Section IIII CI 
the probability po(0 converges to a non-zero value and 
therefore the QW exhibits localisation. 



In the case of odd dimension 2d+ \ we augment the 
coin given by Eq. by the Hadamard coin for the 

extra spatial dimension 

(;j(2d+l) ^(^(2d) ^^(0^Q)_ (58) 

Performing a similar analysis as in the case of even di- 
mensions we find that for the QW driven by the coin 
Q{2d+i) ^Yyq probability that the walk returns to the ori- 
gin decays like due to the Hadamard walk in the extra 
spatial dimension. Hence, this QW is recurrent. 

We note that due the fact that the 2-D Grover walk 
is transient for the initial state ipG the same statement 
holds for the above constructed QWs, supposed that the 
initial state contains tpG in its tensor product decomposi- 
tion. Such vectors form a subspace with dimension equal 
^d-i Jqj, even dimensional walk and 2 x 4*^^^ for odd di- 
mensional walk. 



VI. RECURRENCE OF THE 2-D FOURIER 
WALK 

We now turn to the 2-D Fourier walk driven by the 
coin 




(59) 



As we will see, the Fourier walk does not exhibit local- 
isation. However, the decay of the probability po{t) is 
slowed down to so the Fourier walk is recurrent, ex- 
cept for a subspace of states. 

We start our analysis of the Fourier walk with the ma- 
trix 



UF{ki,k2) - {D{ki) (g> D{k2)) F, 



(60) 



which determines the time evolution in the Fourier pic- 
ture. It seems to be hard to determine the eigenvalues 
of J7i?(fcijfe) analytically. However, wc only need to de- 
termine the saddle points of their phases ujj{ki, k2). For 
this purpose we consider the eigenvalue equation 



$(fci,A:2,w) = det 



(UF{ki,k2) 



0. 



(61) 



This equation gives us the phases U!i{ki,k2) as the solu- 
tions of the implicit function 

$(fci, fc2, = 1 + cos(2fc2) — 2 cos(2a;) 4- 2sin2a;+ 
+4 cos k2 sin uj (sin ki — cos fci ) = . 

(62) 

Using the implicit differentiation we find the derivatives 
of the phase uj 

duj cos k2 sin uj (cos fci -I- sin ki) 

dki cos(2tt') -I- sin(2w) -I- cos k2 cos w (sin ki — cos k2) 

did 2 sin A:2 sinw (cos fci — sin ki) — sin(2fc2) 

dk2 2 (cos(2ti;) -I- sin(2aj) -I- cos k2 cos a; (sin ki — cos fc2)) 

(63) 
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with respect to ki and k2. Though we cannot ehminate 
(jj on the RHS of Eq. we can identify the stationary 



points k 



dki 



0, 



1,2 



(64) 



k=k» 



of w(fci,A:2) with the help of the imphcit function 
$(fci,A:2,w). We find the following: 
{i) a;i^2(fci, ^2) have saddle lines 

71 = (fci,0) and 72 = (fci,7r) 



(zi) all four phases uji{ki, k2) have saddle points for 

fcO^- _^ and kl^±^- 
1 4 4 ^2 



It follows from the case {iii) of Section fill CI that the 
two phases a;i^2(fci, fe) with saddle lines 71.2 are responsi- 
ble for the slow down of the decay of the probability po(i) 
to for the Fourier walk, unless the initial coin state is 
orthogonal to the corresponding eigenvectors ^1,2(^1, ^2) 
at the saddle lines. For such an initial state the prob- 
ability po{t) behaves like t'"^ as the asymptotics of the 
integral given by Eq. ((32|) is determined only by the sad- 
dle points {ii). 

Let us now determine the states V'f which lead to the 
fast decay of the probability that the Fourier walk 
returns to the origin. The states tpp have to be constant 
vectors fulfilling the conditions 



(V'F,t^l,2(k)) =0 Vke7i,2, 



(65) 



which implies that must be a linear combination of 
W3,4(k € 71,2) forming a two-dimensional subspace in Tic- 
For k2 ~ 0,TT we can find the eigenvectors of the matrix 
UF{ki,k2) explicitly 



vi{ki,Q) 


= i'2(fci,7r) 




wi(fci,7r) 


= «2(fcl,0) 




z;3(fci,0) 


= W3(fci,7r) 


= ^(1,0,1,0)- 


■^4(^1,0) 


= Vi(kl,TT) 


= i=(0, 1,0,-1)^ 



The explicit form of Tpp reads 



'ipF{a,b) ^ {a,b,a,~b) , 



(66) 



(67) 



where a, 6 € C. We point out that the particular initial 

state 



, / 1 , 1-A 1 f 1-i 1-i 



-hi 




^= (1,0, CO)'-'" 



X 




40 



80 



(68) 



FIG. 4: Probability distribution after 50 steps and the time 
evolution of the probability po{t) for the Fourier walk with 
the initial state i/j — (1, 0, 0, 0)'^. The upper plot of the prob- 
ability distribution reveals a presence of the central peak. In- 
deed, is not a member of the family ■(//^(o, &). However, in 
contrast to the Grover walk the peak decays vanishes which. 
In the lower plot we illustrate this by showing the probabil- 
ity po{t) multiplied by t to unravel the asymptotic behaviour 



which was identified in [25] as the state which leads to 
a symmetric probability distribution with no peak in 
the neighborhood of the origin belongs to the family de- 
scribed by Eq. ((67| . 

We illustrate the results in Figure [3] and Figure [H] In 
Figure [4] we plot the probability distribution and the 
probability po(^) for the Fourier walk with the initial state 
tp — {1,0,0, 0)^. This vector is not a member of the fam- 
ily tppio-jb) defined by Eq. (p7|) . We find that a central 
peak is present, as depicted on the upper plot. How- 
ever, in contrast to the Grover walk, the peak vanishes 
as shown on the lower plot, where we plot the probabil- 
ity po(^) multiplied by t. Nevertheless, the plot indicates 
that the probability po{t) decays like t^^, in agreement 
with the analytical result. In contrast, for Figure [S] we 
have chosen the initial state given by Eq. ([55]) which is 
a member of the family il)F{o,,b). The upper plot shows 
highly symmetric probability distribution. However, the 
central peak is not present and as the lower plot indicates 
the probability po(i) decays like . 

We conclude that the Fourier walk is recurrent except 
for the two-dimensional subspace of initial states defined 
by Eq. (j67p for which the walk is transient. 
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i / 



/ 



n 




FIG. 5: Probability distribution after 50 steps and the time 
evolution of the probability po{t) for the Fourier walk with the 
initial state given by Eq. (|68p . Since tp is a, member of the 
family ^pp {a,b) the central peak in the probability distribution 
is not present, as depicted on the upper plot. The lower plot 
indicates that the probability po{t) decays like . 



Let us now turn to the estimation of the Polya num- 
bers of the 2-D Fourier walk for the two-dimensional sub- 
space of initial states given by Eq. (p7)) . We make use of 
the normalization condition and the fact that the global 
phase of a state is irrelevant. Hence, we can choose a to 
be non-negative real and b is then given by the relation 



(69) 



Therefore, we parameterize the family of states defined 
by Eq. (|67p by two real parameters — a ranging from 
to and the mutual phase (f) £ [0, 27r). The exact 
expression for po{a, (f>, t) can be written in the form 



PQ{a,cj),t) 



Ki{t) - K2{t)a^\ -a^{cos(f>- 



<2 



(70) 

where ifi.2 has to be determined numerically. Neverthe- 
less, the numerical simulation of po(a,(j),t) at two val- 
ues of (a, 4>) enables us to find the numerical values of 
^1,2 (i) and we can evaluate poia, 4>, t) at any point (a, 4>). 
The probability po(a, t) shows the maximum at a = ^, 



3ir 



phase (p = Consequently, these points also represent 
the maximum and the minimum of the Polya numbers. 




0.314 



n 

4> 



2w 



FIG. 6: Approximation of the Polya numbers for the 2-D 
Fourier walk and the initial states from the family of states 
defined by Eq. (|67|l in their dependence on the parameters 
of the initial state a and (j>. Here we have evaluated the first 
100 terms of po{a, 4>, i) exactly. The Polya numbers cover the 
whole interval between the minimal value of « 0.314 

and the maximal value of p™"^ ~ 0.671. The extreme val- 
ues are attained for a = 1/2 and f^™-'"- = 77r/4, respectively 
^max _ 3.y^^4 Qjj ^jjg lower plot we show the cut at the value 
a = 1/2 containing both the maximum and the minimum. 



In Figure [5] we present the approximation of the Polya 
number Eq. ([5|) in its dependence on a and (j) and a 
cut through the plot at the value a — 1/2. Here we 
have evaluated the first 100 terms of po{a,<j>,t) exactly. 
We see that the values of the Polya number vary from 
0.314 to the maximal value of 
0.671. We note that for the initial states that 
do not belong to the subspace defined by Eq. ([S7)l the 
Polya number equals one. 



the minimum p™"' 



and the minimum for the same value of a and the 



VII. CONCLUSIONS 

Our results, summarized in Table HIl demonstrate that 
there is a remarkable freedom for the value of the Polya 
number in higher dimensions, depending both on the 
initial state and the coin operator, in contrast to the 
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QW 


Section 


Return probability 


Polya number 


1-D unbiased 


IIVAI 


t^^ for any tp 


1 


d-D TP coin 


IIVBI 


for any 'ip 


< 1 

independent of ■)/' 


2-D Grover 


IV A| 


const, for 7^ ^pc 

for ip = ipc 


1 

< 1 


2-D Fourier 


m 


t^^ for ip ^ ipF 

for Tp £ TpF 


1 

< 1 

dependent on t/j 



TABLE II: Summary of the main results. We list the types of 
studied QWs, the asymptotic behaviour of the return prob- 
ability and the Polya number in the respective cases in its 
dependence on the initial state ip. 



classical random walk where the dimension of the lattice 
uniquely defines the recurrence probability. Hence, the 
quantum Polya number is able to indicate physically dif- 
ferent regimes in which a QW can be operated in. We 
expect further interesting effects when we relax the con- 
dition of allowing only for unit steps and introduce larger 
jumps. In that case the size of the coin operator can ex- 
ceed the dimension of the lattice thus in low dimensional 
lattices some effects seen in higher dimensions can be an- 
ticipated, e.g. localisation has been found in three-state 
quantum walks on a 1-D lattice (4l] |. 

In the present paper we assumed a specific measure- 
ment scheme where the dynamics are not continued af- 
ter the measurement is performed. We note that this is 
only one of the possibilities to define the Polya number, 
one could vary the frequency of measurements randomly 
or in a deterministic manner while continuing the time 
evolution. The present definition has the advantage of 
maintaining unitary time evolution, thus a pure state for 
initial pure states. 
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APPENDIX A: RECURRENCE CRITERION FOR 
QUANTUM WALKS 

In this Appendix we show that the recurrence criterion 
for QWs is the same as for RWs, i.e. the Polya number 



equals one if and only if the series 

00 



(Al) 



diverges. 

According to the definition of the Polya number Eq. 
([6]) for QWs we have to prove the equivalence 



-(-00 

t=i 



s = 



(A2) 



We note that the convergence of both the sum S and the 
product P is unaffected if we omit finitely many terms. 

Let us first consider the case when the sequence po (t) 
converges to a non-zero value < a < 1. Obviously, in 
such a case the series S is divergent. Since po {t) converges 
to a we can find for any e > some to such that for all 
t > Iq the inequalities 

1 - a - £ < 1 -po(i) < 1 - a + e- (A3) 
hold. Hence, we can bound the infinite product 

lim (l-a-e)*<P< lim (1 - a + e)* . (A4) 

t — > + oo t — *+oo 

Since we can choose e such that 

|l-a±e|<l, (A5) 

we find that limits both on the left-hand side and the 
right-hand side of Eq. (|A4[) equals zero. Hence, the prod- 
uct P vanishes. 

Let us now turn to the case when poit) converges to 
zero. We denote the partial product 



Pn = iiii-pom 



(A6) 



Since 1 — po{t) > for all ^ > 1 we can consider the 
logarithm 



lnP„ = ^ln(l-poW) 



and rewrite the infinite product as a limit 
P= lim e'"^". 



(A7) 



(AS) 



Since po{t) converges to zero we can find some such 
that for all t > to the value oi po{t) is less or equal than 
1/2. With the help of the inequality 

- 2a; < In (1 - a;) < -a; (A9) 

valid for a; G [0, 1/2] we find the following bounds 



(AlO) 



t=i 



t=i 
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Hence, if the series S is divergent the hmit of the sequence 
(hiP„)^_^ is — oo and according to Eq. (|A8p the product 

P vanishes. If, on the other hand, the series S converges 
the sequence (lnP„)^ is bounded. According to Eq. 

+00 

(|A7|) the partial sums of the series ^ ln(l — po{t)) are 

t=i 

bounded and since it is a series with strictly negative 



terms it converges to some negative value & < 0. Con- 
sequently, the sequence (lnP„)^_^ converges to b and 
according to Eq. (|A8|1 the product equals 

P = e^>0. (All) 

This completes our proof. 
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